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Abstract 

The phase behavior of confined neniatogens is studied using the Lebwohl-Lasher model. For 
three dimensional systems the model is known to exhibit a discontinuous nematic-isotropic phase 



1-^ ' transition, whereas the corresponding two dimensional systems apparently show a continuous 

■ ■ Berezinskii-Kosterlitz-Thouless like transition. In this paper we study the phase transitions of the 

Lebwohl-Lasher model when confined between planar slits of different widths in order to establish 
c/3 \ the behavior of intermediate situations between the pure planar model and the three-dimensional 

C^ ' system, and compare with previous estimates for the critical thickness, i.e. the slit width at which 

the transition switches from continuous to discontinuous. 

O : PACS numbers: 64.60Cn, 61.20.Gy 
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I. INTRODUCTION 



The Lebwohl-Lasher (LL) model l|-|6| is a lattice model of an anisotropic fluid. Each 
site of the lattice is occupied by a uniaxial molecule. A molecule interacts exclusively with 
molecules located at its nearest-neighbor (NN) sites. The total potential energy takes the 
form: 

f/ = -e5^P2(s.-s,), (1) 

where e is the coupling parameter (e > 0), s, and Sj are unit vectors that indicate the 
orientation of the molecules in the corresponding sites, P2 is the second degree Legendre 
polynomial, and < ij > indicates that the sum is restricted to NN pairs of sites. The LL 
model can be deemed as the lattice version of the hard sphere Maier-Saupe (HSMS) fluid 
|7H12| . Most of the simulation work on the LL model has been carried out on simple cubic 



lattices for three-dimensional (3D) systems, and square lattices for the two-dimensional (2D) 
case, although some variations have also been considered |6|. 

A number of papers have been devoted to the analysis of the phase diagram of the LL 
model using computer simulation. The model in 3D has been found to exhibit a discontinu- 
ous nematic-isotropic transition |2|-|5|. The planar Lebwohl-Lasher (PLL) model, defined on 
a square lattice has also been treated extensively using computer simulation 13l-ll7j. From 
this set of results it has been suggested that the PLL model presents a topological defect 
driven continuous transition of the Berezinskii-Kosterlitz-Thouless (BKT) type 18|, [iQJ . No- 
tice however, that some differences between the transition of the PLL model and that of 
the two dimensional XY model (the paradigm for the topological BKT behavior) have been 
recently reported |l7l|. 

In this paper we will pay attention to the nature of the phase transitions of this system 
under confinement in a slit pore, and will study the influence of the pore width on the tran- 
sition. Herein we will be dealing with slits formed by neutral walls, by which the systems 
under consideration will be, in fact, slab models. In this regard, rigorous results 20N22| 
indicate that this type of models cannot support true long range order at finite temper- 
ature (in common with bidimensional systems 231]). This implies that in our context of 
confined/slab and planar systems, we may encounter phases with quasi-long range orienta- 
tional order, which will be here referred to as quasi-nematics. From the point of view of 
simulation, the LL model confined in slit pores was previously studied by Cleaver and Allen 



24j |. They concluded that the system has a critical thickness, He, below which there is no 
bulk-like transition. The existence of such a multicritical point in the T — H plane (where 
T is the temperature and H is the thickness of the slab) can be explained using theoretical 



arguments |25|-|2 Tj . Nonetheless, according to the theoretical approach of Telo da Gania and 
Tarazona|26|. in the case of neutral walls, one should expect He — t- oo. Here, we will address 
this issue resorting to simulation techniques, and analyzing the results obtained for system 



sizes much larger than those considered in Ref . J2j] . We will thus assess the bounds proposed 
therein for such a possible critical thickness. 

In close connection with this work, the effect of the confinement on the isotropic-nematic 
transition has been studied in Ref. 12| for the HSMS model, where it was found that for some 
temperatures the first order isotropic-nematic transition can disappear when the system is 
confined in fiat slits with thickness below a certain width Hc{T). For smaller values of the 
pore width, H < Hc{T), a BKT-like transition appears. Nevertheless, in the HSMS model 
one has to deal with density fluctuations that are not present in the LL model, and this 
could influence the phase behavior. Notice, that it is however possible to introduce density 
fluctuations within a lattice model, as it can be seen in the so-called Lebwohl- Lasher lattice 



gas mo del [28 1 . 



Also related to the present work, a very recent article by Fish and Vink|29|] focuses on 
the effects of conflnement on a generalized version of the LL model, in which the angular 
dependent component of the interaction is oc |sj ■ Sj|^. These authors analyze the behavior 
of the model for values of p > 8 (note that in the present instance p = 2) for which they 
show there is a well deflned critical thickness that vanishes for large values of p when the 
phase transition becomes flrst order even in the two-dimensional limit. 

In summary, when going from the LL bulk behavior to that of the conflned system, 
we should be able to sort out between various possible scenarios. First, the transition 
between that isotropic and nematic phase might be second order, being the ordered phase 
not critical below the transition temperature, with a flnite and non-zero order parameter and 
a diverging susceptibility only at the critical point. This situation is in principle ruled out by 
the exact results that preclude the existence of long range order (i.e. a non vanishing order 



parameter) in our model|20H22|. Another possibility, would be the presence of a continuous 
BKT transition, in which below the transition temperature the system exhibits quasi-long 
range orientational order (quasi-nematic phase with a vanishing order parameter) and the 



susceptibility diverges at all temperatures below the transition temperature. Some subtle 
issues, regarding; what a true BKT transition implies in connection with the discussion of 



Ref. 17| will be addressed in later sections of this paper. Finally, another alternative is 



illustrated by the generalized XY and related models 30|-|34|. which for sufficiently "sharp 



and narrow" interactions J3^ have been shown to undergo a first order transition between 
the isotropic and quasi-nematic phases. It is thus, the aim of this work to provide additional 
information in order to be able to discern between those scenarios. 

The rest of the paper is organized as follows; after this introduction, in section II we 
describe the simulation methodology and summarize the details of the calculations and 
systems under consideration. In section III we present our main results and discuss our 
most relevant conclusions. 

II. SIMULATION TECHNIQUES 

We will deal with systems consisting oi L x L x H sites. Periodic boundary conditions 
(PBC) are applied on the x— and y— directions, and the systems are confined by neutral 
walls in the z— direction. For a given slab thickness, H, results for different values of L are 
taken into account in order to perform the finite-size scaling analysis. We have here studied 
systems with H = 1,2, 3, 4, 5, 8, and 16. For each value of H we have considered a series of 
L values, namely, L = 10, 15, 20, 25, 30, 25, 40, 45, 50, 60, 80, 100, 120, 140, 160 and 200. 

In addition we have also simulated various systems using PBC in the three spatial direc- 
tions. In particular, systems with H = IQ and different values of L, so as to analyze the 
effects of the boundary conditions on the transitions. Fully cubic systems L x L x L with 
PBC were also simulated in order to represent the 3D bulk system. Obviously, we will not be 
dealing here with "true" bulk systems, but we will use the results of non-confined isotropic 
periodic systems, after a finite size scaling analysis is performed, as a good approximation to 
the bulk system results. For simplicity, these corrected results will be referred to as "bulk" 
data. 

We have performed Monte Carlo simulations combining single particle Monte Carlo steps 
with cluster algorithms [4, |9|, |35| using multicluster moves 36|, |33]. For given values of the 



system sizes, L, and H we performed independent simulation runs at several temperatures 
close to the range where the transitions are expected. The results were analyzed using 



efficient re- weighting procedures [36|, l38|. The simulation procedures have been adapted 



from our previous works to the simpler lattice system, and technical details can be found 



elsewhere J9l-lll|. In order to locate the isotropic-quasi-nematic transitions we monitored 



the largest eigenvalue, A_|_ of Saupe's tensor 



1 ^ 

i=l 



(2) 



For a given system size, described by the lengths L and H we can define pseudo-critical 
temperatures, Tc{L,H), in terms of the behavior of A+ as a function of the temperature, 
and also in terms of the temperature dependence of the susceptibility, this quantity being 
defined by means of the fluctuation of the order parameter as: 



X 



N{<Xl>- <X+ >2) /ksT. 



(3) 



In practice, we consider two criteria to determine the pseudocritical points, namely the tem- 
perature at which x, as defined in Eq.Q, is maximum, and the temperature that gives the 
largest value of \dX^/dT\. Then, one can use the pseudo-critical temperatures to extrapo- 
late the transition temperature in the thermodynamic limit {L — t- oo). Follo wing the usual 
practice [l3|, uJ^ , we have used both the expected scaling for BKT transitions 
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43| 



n{L) = n + 
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and the scaling equation of second order transitions 



nL) 



2' 
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l|: 



T,iL) = Te + aiL-i/^ 



(4) 



(5) 



Notice that we have not used the loci of maxima of the excess heat capacity per particle, 
Cy, as an additional criterion to define pseudo-critical temperatures. This alternative was 



used in 



24l ], but most likely is not a good choice for systems that might exhibit a BKT- 



like transition (e.g. for small values of H). In such a case, the maximum in c^ is not well 
defined and does not diverge with increasing sample sizes. Therefore it is not obvious that 
its location signals the presence of a phase transition. It is worth mentioning that we have 
implicitly assumed in Eq. (^ an exponential divergence of the correlation length £_oc e'^^" 



with u = 1/2; this value of u is known to be appropriate for the XY-model|l4j. Il9l . |43 |. 
We decided to use z/ = 1/2 following Ref.(1l4j). due to the fact that a sensible fitting of 
the simulation results to a non-linear equation involving four adjustable parameters would 
require both a much larger range of values of L and very precise simulation results. 
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FIG. 1: (Color on line) Order parameter A-|- as a function of the temperature for different system 
widths, H = 1,4:, 16 and different systems sizes, L. Symbols denote the result of different simulation 
runs, and lines represent the results of the reweighting analysis. The legends in the figures indicate 
the different values of L. 

III. SIMULATION RESULTS 



In Figures [TJIU we depict the temperature dependence of A+, x? dX+/dT and Cy, for 
different system sizes and pore widths. 

It can be seen that the dependence of these properties on L is qualitatively similar for the 
three slit widths considered in the figures. For a given slit width, the susceptibility x diverges 
with L both at the temperature corresponding to the maximum and below. The curves of 
< A_|_ > as a function of T exhibit an inflection point, and the derivative of A_|_ with respect 
to T seems to diverge at a given critical temperature. It is to be stressed that the values 
of A+ below the apparent transition temperature decrease with increasing sample sizes, in 
contrast with the expected behavior from first and second order transitions. On the other 
hand the heat capacity exhibits a maximum which does not diverge with L. The behavior of 
all these properties around the transition temperature resembles that of a topological BKT 
transition, and it is clearly different from what one should expect in the presence of a weak 
first order transition. A second order transition might exhibit non-divergent heat capacity 
curves (with negative a exponent), but the decrease of A+ and divergence of x below the 
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FIG. 2: (Color on line) Order parameter susceptibility, x ^-s a function of the temperature for 
different system widths, H = 1,4, 16 and different systems sizes, L. Symbols and lines as in Fig. [1] 
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FIG. 3: (Color on line) Absolute value of the derivative of the nematic order parameter with 
respect to the temperature, \dX^/dT\ as a function of the temperature for different system widths, 
H = 1,4:, 16 and different systems sizes, L. Symbols and lines as in Fig. [TJ 

pseudo-critical Tc{L, H), fit better into the picture of a continuous phase change which shares 
a number of features with the continuous BKT transition. For the sake of comparison, in 
Figure [5] we summarize the results of simulations for unconfined systems using cubic boxes of 
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FIG. 4: (Color on line) Reduced excess heat capacity as a function of the temperature for different 
system widths, H = 1,4:, 16 and different systems sizes, L. Symbols and lines as in Fig. [H 
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FIG. 5: (Color on line) Temperature and system size dependence of various properties of the 3D LL 
model (cubic cells without walls). The length of the cell side, L, is shown in the legends. Symbols 
and lines as in previous figures. 

different sizes with full PBC. It seems evident that the qualitative behavior of the confined 
system is quite different from that of the bulk, which is known to present a weak first order 



isotropic-nematic transition 



i- 



Returning to the confined system, in table [T] we gather the results for the estimates of its 
transition temperatures for different slit widths calculated using the two aforementioned def- 



initions of the pseudo-critical temperatures, and the two scahng laws. The results for H = 1 



agree with those reported in Ref . [ij] , but differ slightly from those reported in Ref . 13| using 
the scaling laws of second order transitions. The results of the table show that the estimates 
of the transition temperature are conditioned by the scaling law used in the extrapolation 
to the thermodynamic limit. However, the results hardly depend, within error bars, on 
the particular definition of the pseudo-critical temperature. The variation of the transition 
temperature with H is monotonic, and the transition temperatures approach smoothly the 
bulk value as H increases. This is more clearly seen in Figure [6l where one can appreciate 
the quasi- linear dependence of Tc (as calculated from (^) on 1/H. This Kelvin-like scaling 
of the transition temperature leads to an extrapolated value lim//_j.oo Tc{H) = 1.123 ± 0.005 
that agrees rather well with the bulk value T^ = 1.1225(1), which we have obtained using 
cubic systems with full PBC, and in accordance with the results of Priezjev and Pelcovits|4]. 
Note, however that the Kelvin scaling only applies strictly to first order phase transitions. In 



the case of second order transitions, correction terms must be incorporated 40l-l42| . From 
our discussion it is clear that in our case a first order phase transition is ruled out, so 
deviations from linearity could in principle be ascribed to the continuous character of the 
transition. It is worth stressing that T^ estimates become independent of the scaling relation 
used as H increases. This is an indication, that even if the transition can still be cast into 
the BKT-like type for growing H, its scaling behavior is gradually switching to that of a 
regular order- disorder transition. 

We also include in Table [T] the estimates of the scaling exponent for the maximum of the 
susceptibility, 7/1^, which can be drawn from the scaling relation: 

X'^''%L) ~ r^". (6) 

The results of this exponent depend on the pore width, and constitute further evidence that 
no first-order transition appears for the system sizes considered. One should expect in this 
latter instance a scaling of the type x™'"^(L) ~ HL^, well away from the values obtained 
here for any width. Incidentally, in the H = 1 case the value is relatively close to the 



two dimensional Ising critical exponent |44j|. 7/z/ = 7/4, and in agreement with the value 
reported by Kunz and Zumbach 1^ 7/z/ = 1.72 ± 0.05. For larger slit widths, the value of 
7/z/ decreases, what further deviates from the limiting behavior of a first order transition 
when // — 7- 00 (7/z/ = 2). This implies that in the range < if < 00 one should expect a 
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FIG. 6: (Color on line) Pore width dependence of the transition temperatures estimated using the 



scahng law (g]). T 
Refs.y, y, and 



re result of a linear fit to 1/H is represented by a solid line. Values taken from 



li 



are also included for comparison. Note that the bibliographic values and 



those of this work fall on top of each other and can hardly be distinguished. 



121 and 



non-inonotonic behavior of 7/z/, as was already found in the confined HSMS fluid 
it is a clear indication that if = 16 is still far away from the first order transition limit. 
This situation is in contrast with the results recently reported by Fish and Vink[29j for the 
generalized LL model with p = 8. In this case the angular interaction is much narrower than 

32|,|34| 



that of the simple LL model and bears some resemblance with the g-state Potts model 
(with q oc ^/p)■ Fish and Vink found that '-//u grows from 1.63 in the two dimensional limit 
approaching 7/// — )■ 2 as the critical He is reached and the continuous transition develops 
into a first order transition. We assume that as p decreases He increases (as observed in 
Ref. [29| for p > 8), to the point that for p = 2 the determination of the critical thickness is 
well beyond our present computational capabilities. On the other hand, it is to be noticed 
the fairly regular dependence of Tc{H) on the pore with. 

From all this evidence, and in particular, from the size dependence of the order parameter 
and the susceptibility, one can conclude that the isotropic-quasi-nematic transitions found 
for all the confined systems under scrutiny {1 < H < 16) do not fit in the picture of 



10 



TABLE I: Estimates of the isotropic-quasi-nematic transition temperature in the thermodynamic 
hmit, using different prescription of the pseudo-critical temperatures and scahng laws (See the 
text for details); and scaling exponents 7/1^ for the maxima of the susceptibility x- The transition 
temperature in the bulk system is Tn-i = 1.1225(1). Error bars across the table are shown between 
parentheses in units of the last figure and correspond to a confidence level of 95 % 



H 


1 2 3 4 5 8 16 


Tc [X, Eq.©] 


0.514(2) 0.772(4) 0.889(3) 0.963(3) 1.000(2) 1.062(1) 1.104(1) 


Tc [\dX+/dT\, Eq.dH)] 


0.508(4) 0.764(7) 0.889(5) 0.956(4) 0.999(3) 1.062(2) 1.104(1) 


Tc [X, Eq.®] 


0.536(3) 0.782(3) 0.905(3) 0.973(3) 1.010(2) 1.067(1) 1.105(1) 


T, [\dX+/dT\,Eq.m 


0.531(4) 0.786(5) 0.906(4) 0.969(3) 1.008(2) 1.067(1) 1.105(1) 


-y/u 


1.69(2) 1.65(3) 1.63(3) 1.59(3) 1.60(4) 1.49(3) 1.40(7) 



first or second order phase changes, and share some resemblance with the continuous BKT 
transition. The results also indicate that the possible critical thickness He of the Lebwohl- 
Lasher model, if exists, must appear for He > 16. Therefore, the value of He reported 
in Ref. [24] is most likely underestimated. It is possible to further analyze the effects 
of dimensionality on the transition if the walls in the z— direction are replaced by PBC, 
but still dealing with the 2;-direction on a different footing as compared with the x— and 
y— directions. More precisely, we will consider a series of systems oi L x L x H sites 
with PBC on all three directions and for H = 2,4, and 16. This anisotropic LL model 
is expected to enhance the correlations of the nematogen orientations in the 2;— direction, 
and eventually lead to the phase behavior of the bulk system as predicted by theoretical 
arguments 25|. Interestingly, we have found that the anisotropic LL model with H = 2, and 
H = 4 also exhibits BKT-like transitions similar to those of the corresponding confined LL 
system occurring at slightly higher temperatures. Moreover, the same behavior is found for 
if = 16; the system with PBC clearly shows a dependence of c^"^ with L inconsistent with 
a first order transition. This can be appreciated in Figure [TJ where we present the results for 
the value of the maximum of the excess heat capacity per molecule, c™"^ for both systems. 
It can be seen that for both confined, and PBC systems c^"^ does not diverge. In the same 
Figure we include the result for cubic systems H = L with PBC (bulk LL model); in this 
case the expected scaling behavior, d^""^ ~ A^, of a first-order transition is observed. From 
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FIG. 7: (Color on line) Value of the maximum of excess heat capacity per particle as a function 
of the system size for cubic systems (H=L) with PBC (bulk LL), systems with H=16 and PBC 
(anisotropic LL), and systems with H=16 confined between neutral walls. The line connecting 
points for the bulk system corresponds to a fit of the results to the equation Cv{L) = ao + aiL^ 



the results of c^"-^ it is possible to compute the latent heat, AE of the transition 45| : 



^<-)--i|.-^ 



(7) 



where c± is related with the specific heats of the two phases^]. By fitting the results for 
L > 25 we get AE/e = 0.0584 ± 0.0013. 

In the discussion above, we have purposefully avoided an explicit reference to the ques- 



tions recently raised by Paredes et al 17|, |46| as to the existence of a true BKT transition 
in the PLL model. After performing a finite-size scaling analysis of the simulation results 



at temperatures around anc. 
literature, Paredes et al 17, 



be 



ow the estimates of the transition temperature found in the 



461] conclude that the PLL lacks a true topological transition. 



Using results for several of system sizes, they infer that the L-dependence of the order pa- 
rameter distribution for T < Tbkt does not follow the expected scaling for a line of critical 
points. In particular, they argue that the lack of crossing of the Binder cumulant [45 
curves for different system sizes at a fixed temperature is a strong evidence of the absence 
of quasi long range order in the PLL model. In order to gain some extra insight into this 



47 



48| 



12 



problem, we have performed series of simulations with a broad range of system sizes at three 
temperatures: T* = 0.50 (slightly below the range of our Tc estimates), T* = 0.54 (slightly 
above), and T* = 0.60. The so-called Binder cumulant, g^ can be defined as(48|: 



< ^t > 

^4 = ^TT^- (8) 



<A1>2- 



It is well established 



45j that for second order transitions at the critical temperature g^ 
reaches (for large system sizes) a critical value (different from those corresponding to ordered 
and disordered phases) which becomes independent on the system size, gi{L,T^ = g^^' . 
According to the usual description of the topological transitions, below Tbkt there should 
be a line of critical points, and therefore at a fixed temperature g/^{L\T) should approach 
a critical value gl{T) as L — )■ oo. This value must be different from those of the isotropic 
{g4 ^ 3/2) and nematic ((74 ~ 1) phases. From this point of view one can expect that plotting 
g4 as a function of T, the curves with different values of L should merge for T < Tbkt- Of 
course, finite size effects could eventually lead to a small degree of crossing (See Ref. (iTj). 
Therefore, from our point of view the absence of crossing between the gi{T) curves with 
different values of L must not be regarded as a signature of lack of criticality. In figures [8] 
and [9] we show the results of g^ as a function of the system size for T* = 0.50 and T* = 0.54, 
and T = 0.60. The results seem to be compatible with the presence of a BKT-like transition 
in the PLL model. For large system sizes, at T = 0.54, and T = 0.60 g^iL) seems to 
approach to the expected value for the isotropic phase (5^4 ~ 3/2), whereas for T = 0.50 
(with system sizes up to L = 896) the values of g^i^L) apparently converge towards a critical 
value as L — 7- 00. 



Another point raised by Paredes et al. 17|] concerns the apparent violation of the hyper- 



scaling relation of the critical exponents, 2/3/i/ + 7/z/ = d (where d is the space dimensional- 
ity). In the case of a BKT transition, the exponent z/ is not defined, but the exponent ratios 



can still be calculated 



491]. According to Ref. 17|] in the case of the PLL this relation is only 



fulfilled within a 3% accuracy, one order of magnitude less than in the case of the XY-model. 
In our case, calculations carried out at T* = 0.50 (below the transition temperature) also 
indicate deviations around 5%, somewhat larger than the statistical uncertainties. Inter- 
estingly, previous calculations performed at the transition temperature for the continuum 



HSMS model Uj agree with the hyperscaling behavior within a 0.7% error. Moreover, using 
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FIG. 8: (Color on line) Cuniulant ratio 54 =< A^ > / < A^ >^ as a function of the inverse of the 
system size for T* = 0.50, 0.54, and 0.50 
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FIG. 9: (Color on line) Detail of the Cumulant ratio g'4 =< A^ > / < A^ >^ as a function of the 
inverse of the system size for T = 0.50 



an alternative definition of the susceptibility for temperatures above Tc{H, L)[45|, l50|, l51| 



X = N<\1> /ksT, 

we found that the hyperscaling relation is appropriately fulfilled. 

Some additional information can be obtained from an analysis of the percolation of the 
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FIG. 10: (Color on line) Probability of finding cluster percolation in, at least, one direction in 
the simulation of the three-dimensional Lebwohl-Lasher model as a function of the temperature. 
Different curves represent results for different system sizes. 



clusters constructed by the simulation algorithm as a function of the temperature, so as 
to evaluate the degree of correlation between the particle orientations within the simulated 
samples. Let us recall that the Swendsen- Wang-like (SW) algorithm applied in this work 
belong to the class of rejection-free cluster methods, and for some simple systems the tem- 
perature at which the cluster percolation occurs corresponds to that of the phase transition. 
This property was used by Tomita and Okabe 43|] to locate the BKT transitions of two- 
dimensional XY and Potts-Clock models in two dimensions. For the PLL model we have 
carried out multi-temperature simulations using the single tempering algorithm of Zhang 
and Ma[52|. In Figure [TO] we present the results of the percolation probability, Xper, defined 
as the fraction of configurations containing at least one percolating cluster, for the 3D LL 
model with PBC. It can be seen that the percolation threshold appears at a temperature 
slightly above the nematic-isotropic transition temperature. In addition. Xper{L, T) shows a 
non-monotonic behavior with T for large system sizes at temperatures close to the thermo- 
dynamic transition. The behavior of Xper{L,T) is qualitatively similar for the PLL model 
(See Fig. [TTj) : and at temperatures close to the Tbkt estimates the curves for different 
system sizes show a clear tendency to merge. The crossing of the curves for different system 
sizes observed for large system size seem to indicate that the aforementioned merging is not 
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just a consequence of correlations induced by the periodic boundary conditions. Moreover, 
in a similar percolation analysis carried out by us for the planar Heisenberg model (3D spins 
on a plane), the Xper{L,T) curves for different system sizes do not exhibit any crossing at 
finite temperatures, but seemingly merge as T — t- 0. This supports the general view 33| that 
the 2D Heisenberg model does not have a phase transition at T > 0, and underlines the 
essentially different phase behavior of the PLL model. In the thermodynamic limit. Figure 
[To] seems to indicate that in the case of the 3D LL model, there should be an abrupt change 
from the non-percolating state {Xper = 0) to a fully percolating state {Xper = 1) at a finite 
temperature, which fits into the picture of a first order transition between the isotropic and 
a truly nematic phase. In contrast, in Figure [HI one finds that in the PLL model, at least 
for the system sizes here considered, the fully percolating state is only reached at T = 0. 
Note, that by construction, the SW cluster algorithm may yield Xpgj. < 1 for finite tempera- 
tures even in the case of truly orient ationaly ordered states. The presence of the maximum 
after the first crossing (occurring both in the 3D LL and PLL models for finite sizes, but 
seemingly not present in the 3D Heisenberg system) might then well be an effect of the 
cluster algorithm. On the other hand, analysing the size dependence of the curves plotted 
in Figure [Til one is tempted to assume that the maxima will continue to grow and shift to 
lower T as L increases, until finally Xper = 1 is reached for a given sample size. Whether 
this is really the case, and if so, Xp^r = 1 is reached at T > or not, cannot be assessed at 
present using reasonable computer resources. 

In any case, we believe that the percolation analysis sketched above confirms that the 
PLL model indeed presents a phase transition. It might be the case, that we are not dealing 
here with a strict BKT transition, if one takes into account the previous discussion on the 
hyperscaling relation, but its phenomenology is closely related to that of the BKT transition. 
On the other hand, the other anomalies in the model's scaling behavior found by Paredes 
et al. [17| could be ascribed to finite size effects. To conclude this discussion about the 
likelihood of a topological transition for the PLL model, it is worth to compare the phase 
behavior of the three-dimensional XY and LL models. The three-dimensional XY model 



presents a continuous transition 53| without a divergence in the specific heat, whereas the 



three-dimensional LL model exhibits a first order transition. Taking as a reference the 



critical behavior of Potts models 



we would not expect in principle that the PLL model 



had a weaker transition than that of the two-dimensional XY model. 
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FIG. 11: (Color on line) Probability of finding cluster percolation in, at least, one direction, in the 
simulation of the planar Lebwohl-Lasher model as a function of the temperature. Different curves 
represent results for different system sizes. 

In summary, we have studied the order- disorder transition for the confined LL model by 
means of Monte Carlo simulation and finite-size scaling analysis. Our results indicate that 
the critical pore width signaling the crossover between bulk and 2D behavior must be larger 
than the values indicated by previous simulations. The need for a reliable finite-size scaling 
analysis on systems with larger widths, which would imply simulations for much larger 
systems hampers the estimation of He. In addition, our results for slit-like systems with 
full PBC (anisotropic LL model) suggest that if H^ has a finite value, it will likely be much 
larger than H = 16. Moreover, the fact that the critical exponent relation 7/z/ for the pore 
widths considered does not yet show any trend to converge towards the expected behavior in 
a first order transition, is a further indication that we are very likely away from the critical 
thickness. This might fit into the picture drawn by Telo da Gama and Tarazona 26l| . who 
suggest that one should expect He ^ 00. However, in Ref. l26| it is argued that spin- waves 
would destroy the ordered phase for any finite H, but no BKT transition would occur. Our 
findings suggest that an order-disorder phase transition with some BKT-like features does 
indeed take place, as the divergence of the susceptibility, its size dependence, the crossover 
of percolation curves and the size dependence of the order parameter seem to evidence. It 
is worth pointing out that the situation depicted here is in marked contrast with the abrupt 
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switch from continuous 2D melting behavior to discontinuous first order melting which has 
been ar gue d to occur in hard sphere colloidal models when going from monolayer to bilayer 



systems 



55|. 
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